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ABSTRACT 

Context. The secondary-to-primary boron-to-carbon (B/C) ratio is widely used to study the cosmic-ray (CR) propagation processes 

in the Galaxy. It is usually assumed that secondary nuclei such as Li-Be-B are generated entirely by collisions of heavier CR nuclei 

with the interstellar medium (ISM). 

Aims. We study the CR propagation under a scenario where secondary nuclei can also be produced or accelerated by Galactic sources. 

We consider the processes of hadronic interactions inside supernova remnants (SNRs) and the re-acceleration of background CRs in 

strong shocks. We investigate their impact in the propagation parameter determination within present and future data. 

Methods. Analytical calculations are performed in the frameworks of the diffusive shock acceleration theory and the diffusive halo 

model of CR transport. Statistical analyses are performed to deteiTnine the propagation parameters and their uncertainty bounds using 

existing data on the B/G ratio, as well as the simulated data expected from the AMS-02 experiment. 

Results. The spectra of Li-Be-B nuclei emitted from SNRs are harder than those due to CR collisions with the ISM. The secondary-to- 

primaiy ratios flatten significantly at ~ TeV/n energies, both from spallation and re-acceleration in the sources. The two mechanisms 

are complementary to each other and depend on the properties of the local ISM around the expanding remnants. The secondary 

production in SNRs is significant for dense background media, «i > 1 cm"', while the amount of re-accelerated CRs is relevant to 

SNRs expanding into rarefied media, rii < 0.1 cm"'. Owing these effects, the diffusion parameter S may be underestimated by a factor 

of ~ 5-15%. Our estimations indicate that an experiment of the AMS-02 caliber can constrain the key propagation parameters, while 

breaking the source-transport degeneracy for a wide class of B/C-consistent models. 

Conclusions. Given the precision of the data expected from ongoing experiments, the SNR production/acceleration of secondary 

nuclei should be considered, if any, to prevent a possible mis-determination of the CR transport parameters. 

Key words, cosmic rays — supernova remnants — acceleration of particles — nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

The problems of the origin and propagation of the charged cos- 
mic rays (CRs) in the Galaxy are among the major topics of re- 
sarch in modern astrophysics. It is generally accepted that pri- 
mary CR nuclei such as H, He, C, N, and O, are accelerated 
in supernova remnants (SNRs) via diffusive shock acceleration 
(DS A) mechanisms, that produce power-law momentum spectra 
(Drury, 1983). At relativistic energies, S oc p^^ ~ E^^ . After 
being accelerated, CRs are released in the circumstellar environ- 
ment, where they diffuse through the turbulent magnetic fields 
and interact with interstellar matter (ISM) (Strong et al., 2007). 
Owing to diffusion, CRs stream out from the Galaxy on a char- 
acteristic timescale Tesc °^ E^''. The spectrum of primary CR 
nuclei predicted at Earth is therefore Np ~ 5 Tesc- The collisions 
of these nuclei with interstellar gas are believed to be the mecha- 
nism producing the secondary CR nuclei, such as Li, Be, and B, 
which are under-abundant in the thermal ISM. Thus, the equilib- 
rium spectra of secondary CRs are E^^ times softer than those 
of their progenitors. At energies above some tenths of GeV per 
nucleon, where CR nuclei reach the pure diffusive regime, this 
picture predicts power-law distributions such as Np ~ E^^^'' for 
primary nuclei and A^s/A^y, ~ £"* for secondary-to-primary ratios 
at Earth. These trends may be straightforwardly derived from 
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the analytical solutions of Maurin et al. (2001). Present obser- 
vations indicate that S ~ 0.3-0.7 and v ~ 2.0-2.4. The bulk of 
the data is collected at E < 10 GeV nucleon"', where the CR 
spectra are shaped by additional effects such as diffusive reac- 
celeration, galactic wind convection, energy losses, and solar 
modulation. Since there is no firm theoretical prediction of the 
key parameters associated with these effects, it is very difficult 
to distinguish each physical component using the experimental 
data. The boron-to-carbon (B/C) ratio is the most robustly mea- 
sured secondary-to-primary ratio and is used to constrain several 
model parameters. Throughout this paper, we call secondaries 
all CR nuclei produced by hadronic interactions, independently 
on the place of origin. The standard approach, hereafter refer- 
ence model, assumes that the secondary nuclei are absent from 
the CR sources. 

In this paper, we examine two mechanisms producing a 
source component of secondary CRs: (i) the fragmentation 
of CR nuclei inside SNRs and (ii) the re-acceleration by 
SNRs of pre-existing CR particles. The secondary CR pro- 
duction inside SNRs was studied in Berezhko et al. (2003) 
and subsequently reconsidered to describe the positron frac- 
tion (Blasi, 2009; Ahlers et al., 2009). Predictions of the p/p 
ratio (Blasi & Serpico, 2009; Fujita et al., 2009) and B/C ra- 
tio (Mertsch & Sarkar, 2009; Thoudam & Horandel, 201 1) have 
also been investigated. An interesting aspect of this mechanism 
is that, if the secondary fragments start the DS A, the secondary- 
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to-primary ratios must eventually increase. Similarly, the re- 
acceleration of background CRs interacting with the expand- 
ing SNR shells may induce a significant transformation of their 
spectra at high energies (Berezhkoet al., 2003; Wandeletal., 
1987). In particular, the re-acceleration redistributes the spec- 
trum of secondary nuclei to a spectrum S ~ £ ''. The main 
feature of both mechanisms is that they produce harder spectra 
of secondary nuclei than in the case of their standard production 
from primary CR collisions in the ISM. These source compo- 
nents of secondary CRs may become relevant at ~ TeV energies. 
Thus, disregarding these eff'ects may lead to a mis-determination 
of the CR transport parameters. The aim of this paper is to ex- 
amine their impact on the CR propagation physics. This task re- 
quires a description of the CR acceleration processes in SNRs 
and their interstellar propagation. In this work, we use fully an- 
alytical calculations in the frameworks of the linear DS A theory 
and the diffusion halo model (DHM) of CR transport. In Sect. §2, 
we present the DSA calculations for CR nuclei, including stan- 
dard injection from the thermal ISM, hadronic interactions, and 
re-acceleration. Sect. §3 outlines the basic elements of the DHM 
galactic propagation. In Sect. §4, we show our model predictions 
for the CR spectra and ratios at Earth, and study the impact of 
the secondary source components on the determination of the 
CR transport parameters. In Sect. §5, we present our estimates 
for the Alpha Magnetic Spectrometer (AMS). We present our 
conclusions in Sect. §6. 



2. Acceleration in SNR shock waves 

We compute the spectrum of CR ions accelerated in SNRs us- 
ing the DSA theory (Drury, 1983), including the loss and source 
terms and both the production and acceleration of secondary 
fragments. Our derivation is formally similar to that in Morlino 
(2011), but the physical problem is the same to that treated in 
Mertsch & Sarkar (2009). Within this formalism, we also com- 
pute the re-acceleration of pre-existing CR particles. 



2.1. DSA Calculations 

We consider the case of plane shock geometry and a test-particle 
approximation, i.e., we ignore the feedback of the CR pressure 
on the shock dynamics. The shock front is in its rest-frame at 
X = 0. The un-shocked upstream plasma flows in from x < at 
speed Ml (density «i) and the shocked downstream plasma flows 
out to X > at speed U2 (density «2)- These quantities are related 
by the compression ratio r - uijuj - na/ni- The particle spectra 
are described by the phase space density f{p, x). The equation 
that describes the diffusive transport and convection at the shock 
for a y'-type nucleus (charge Zy and mass number Aj) is given 
by: 



dx ^' dx^ i,dx^ dp ' ^' ^'' 



(1) 



where Dj{p) is the diffusion coefficient near the SNR shock, u 
is the fluid velocity, P"''' = /Sjcncr^'^^^ is the total destruction rate 
for fragmentation (see Sect. §2.3), cr'"''^ is the cross-section for 
the process, and Qjix, p) represents the source term. Solutions of 
Eq. I can be found, separately, in the regions upstream {x < 0) 
and downstream (x > 0) of the shock front, by requiring that 
df/dx - for X — > +oo. We drop the label j characterizing the 
nuclear species, and make use of the subscript / = I (; = 2) to 
indicate the quantities in the upstream (downstream) region. We 



define the quantities 



D 



(A; - 1) 



Di 



r(A/ + l). 



(2) 



where A, - Jl + 4D,T™V"^ The solution can be expressed in 
the form 



(Mp)e 



fix,p) = 



"lAi 



/o(p)e 



+ 4-12-V , 



U2 + V2 + W2 



(x<0), 



{x>0), 



(3) 



where the downstream integral terms U2, V2, and W2 are given 
by: 

X+oo 
Q2(x' , p)e^-''-^"-'yx' , 



U2(X. 

V2C 



x,p) = + \ Q2{x',p)e-'^^^^'-^'^dx' , 
Jo 

P+00 

Q2(x\p)e-^-^^'-'^''-''^dx' . 
Jo 



(4) 



W2{x,p) 



In the upstream region, U\, Vi, and Wi are still given by Eq. 4 
after performing the substitutions 1 — > 2, /C2 ^ —M, ^2 -^ -ki, 
and 00 — » -00. The distribution function at the shock position, /o, 
is determined by the matching conditions at x = 0. We integrate 
Eq. 1 in a thin region across the shock front. Assuming that D = 
D\ = D2, we find the equation for /o 



P^- = -afoip) - ajip) + —G{p) , 
op u\ 



(5) 



where a — 3ui/(ui - U2) is the known DSA spectral index. The 
term G denotes the sum of the upstream and downstream source 
integrals 



G{p) 



f 

•J —Cx 



gie^'ii-^Vx' + I 026-^"'-'' dx' . 



(6) 



The function j{p) is linked to the destruction term F""''. It is 
defined as 

7(/')=^(Ai-l) + :^(A2-l). (7) 

2 2r 



(8) 



After defining the function 

C jip") 

X(p,p')^a — —'^P"^ 
Jp' P 

the solution of Eq. 5 can be expressed in the simple form 

/o(p) = ar(4^^.-*-'>^. (9) 

Jo \P I Ml P 

From Eq. 3, one recovers the standard DSA solution by setting 
pinei _ Q ^jjQ interactions) and assuming that the injection oc- 
curs only at the shock front (t/, - V, = W, = 0): one finds 
that /2 = /() and /i - foe"'^^'^', while Eq. 9 gives a spectrum 
p ", provided that the source term G(p) is softer than p " (see 
Sect. §2.4). Some simplifications can be made by analyzing the 
timescales of the problem. The DSA acceleration rate for par- 
ticles of momentum p in a stationary shock is F"'^'^ ~ ^^, 
where we assumed Bohm diffusion {D oc p) and strong shocks 
(r a 4). For a SNR of age Tsm-, the condition F^'^'^ s t^^^. de- 
fines the maximum momentum p™™ attainable by DSA. In the 
presence of hadronic interactions, the requirement F'™' «: F"'^'^ 
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must be fulfilled. These relations imply that 20 P'^'^'D/m^ «c 1 
and xP"'^V" ^ 1 at all the energies considered. Under these 
conditions, we can linearly expand A a; 1 + 2P™'D/m-, so that 
A «j 2P°'=V« and k ^ lujD + 2P"'=Vm- The exponential terms of 
Eq. 3 and Eq. 4 can also be expanded as e^'^ ^ \ — ^—x and 
g2« ~ e"^^ (1 + ^xj. Thus, the function j{p) of Eq. 7 is given 
by 



Up) ^a{\ + r) 

and the integral of Eq. 8 by 



P-'D(p) 



""■inel 



X{p,p') ^ a{l + r')^[D(p) - Dip')] . 



(10) 



(11) 



which recovers the expression of Mertsch & Sarkar (2009). 
Below we present the DSA solutions for primary nuclei (in- 
jected at the shock), their secondary fragments (generated in the 
SNR environment), and pre-existing CR particles that undergo 
re-acceleration. 



2.2. Acceleration of primary nuclei 

The injection of ambient particles is assumed to occur imme- 
diately upstream from the shock at momentum /?'"J. The source 
term for primary nuclei is 



QP''(x,p)^Y6(x)6(p-p'">). 



(12) 



Particles can be injected only when their Lannor radius is large 
enough to cross the shock thickness. Thus, we assume a refer- 
ence injection rigidity for all nuclei, R^'^\ so that p'^J = Z/?'"J. 
The constant Y reflects the particle abundances in the ISM and 
their injection efficiencies. In this work, they are determined by 
the data. The phase space density profile is given by 



[Mp)e"^'l'^(P^ (x<0), 

yhip) (^1 - -S7JCJ {x>0). 



(13) 



The upstream profile, ~ e"' '^^, indicates that the plasma is con- 
fined near the shock within a typical distance ~ Dlu\. Owing to 
advection, the particles are accumulated in the downstream re- 

gion, where destruction processes give rise to the term ^— x. The 
momentum spectrum at the shock position is /o oc e^^p^", which 
is the known DSA power-law behavior times an exponential fac 
tor, ~ e-^, given by;r ~ cP'^^'/r"". The condition r"'='= 
implies that;^' < 1. 



« Tsnr 



is implicitly summed over the circumstellar abundances (the hy- 
drogen and helium components). We have assumed that the ki- 
netic energy per nucleon is conserved in the process, i.e., the 
fragments are ejected with momentum ^,ij = Aj/Ak times smaller 
than that of their parents. The presence of short-lived isotopes 
{ghost nuclei) such as 'Li or ''C is reabsorbed in the definition 
of (Titj, while long-lived isotopes such as '°Be or ^^Al (lifetime 
~ 1 Myr) are considered as stable during the acceleration pro- 
cess (timescale Tsnr ~ lOkyr). The source spatial profile takes 
the form of the progenitor nucleus of Eq. 13. In the upstream 
region, one has Q "^ - qi^kje"''^^'^'', where D^ is the diffusion co- 
efficient of the progenitor nucleus at the momentum pl^kj- For 
D oc /7/Z, one can write Dkipl^kj) = (kjDjip), where (kj = j^.. 
In practice, ^^j ~ 1 for all processes k ^> j with Zy,A > 2. The 
terms at the shock, qy^t and qijk, are given by 



^-k]hk{pi^kj)f[;tj^ 



qi,kj{p) 

and the downstream solution reads 

qi.jip) 



fi.jix, p) = fojip) 



U2 



"TfhjiP) 



X, 



(15) 



(16) 



where foj is given from Eq. 9 using G{p) by the expression 

Gkjip) = qukjip)^^ Ukl + r^) . (17) 

We solve all equations starting from the heaviest element and 
proceeding downward in mass. We are interested in the total 
contribution of SNRs to the Galactic CR population, which we 
evaluate to be the integral of the downstream solution over the 
SNR volume left behind the shock 



?snr r 

Jo 



Sj\p) ^ Anp"R,,, 47:x'f2j(x,p)dx, 



(18) 



where x,nax - «2Tsm and '^Jsni is the supernovae explosion rate 
per unit volume in the Galaxy. 

2.4. Re-acceleration of bacl<ground CR nuclei 

Together with the thermal ISM particles of Sect. §2.2, SNR 
shock waves may also accelerate the background CRs at equi- 
librium (Berezhko et al., 2003; Wandel et al., 1987). We refer to 
this mechanism as the re-acceleration of CRs in SNRs. For a 
prescribed distribution function of background CRs, / "(p), the 
DSA solution at the shock is simply 

foM)-^r('-\frip')'4- (19) 



2.3. Production and acceleration of secondary nuclei 

Secondary nuclei originate in the SNR environment from the 
spallation of heavier nuclei on the background medium. The 
source term for a y'-type CR species arises from the sum of all 
heavier k-type nuclei, Q'^^'^ = Tjk>i Qir^- Each partial contribu- 



tion is defined to be: 



k>j '^kj 



1 r°° 

Q'k]Hx,p) = ^^T J ,„. Nk(x,p')f^]Hp')S{p - ^kjP')dp' , (14) 

where Nk(x,p') - A-np' fk(x,p') is the progenitor number den- 
sity and r^'^^ - /^kcncr^''^ is the k ^ j fragmentation rate, which 



Assuming, for illustrative purposes, a power-law form f-^{p) - 
Yj{plp^fy\ the resulting re-acceleration spectrum is 

Since the CR equilibrium spectrum, f^^ oc p^\ is softer than the 



test-particle one {s > a), for p :s^ p . one obtains 



f^fp^-^^{pipn 



(21) 



That is, the effect of re-acceleration is to re-distribute the CR 
spectrum to p^" . Interestingly, in the opposite case {s < a) 
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Fig. 1. Energy spectra of the CR elements Li, Be, B, C, N, O, Ne, Mg and Si. The solid lines represent the reference model prediction. The 
dashed lines indicate the secondary CR component arising by collisions of heavier nuclei in the ISM. The model parameters are listed in Table 1. 
Data are from HEA03-C2 (Engelmann et al., 1990), CREAM (Ahn et al., 2009), AlVIS-Ol (Aguilar et al., 2010), TRACER (Ave et al., 2008; 
Obermeier et al., 2011), ATIC-2 (Panov et al., 2009), CRN (Muller et al, 1991), Simon et al. (1980), Lezniak & Webber (1978) and Orth et al. 
(1978). The Li-Be-B data from CREAM and AlVIS-Ol are combined with our model to obtain the spectra from their secondary-to-primary ratios. 



the re-accelerated spectrum maintains its spectral shape p ■', 
while its normalization is amplified by the factor al{s - a). 
In our model, however, the background spectrum f'^ is com- 
puted as discussed in Sect. §3 and takes the SNR spectra as in- 
put. Therefore, Eq. 19 is an integro-differential equation where 
the DSA-mechanism is fed by its DHM-propagated solution and 
vice-versa. On the other hand, the bulk of the re-accelerated CRs 
come from the low-energy part of the spectrum (below ~ 10 GV 
of rigidity), where the equilibrium CR spectra are fixed by the 
observations so they cannot vary too much. It can be safely as- 
sumed that re-accelerated CRs are a sub-dominant component 
of the total (integral) flux. Hence, we proceed using an iterative 
method as outlined in Sect. §4.5. 

3. Interstellar propagation 

We use a DHM to describe the CR transport and interactions in 
the ISM in a two-dimensional geometry. We disregard the ef- 
fects of energy losses, diffusive reacceleration, and convection. 
The Galaxy is modeled as a disk of half-thickness h, containing 
the gas and the CR sources. The disk is surrounded by a cylin- 
drical diffusive halo of half-thickness L, radius rm^x, and zero 
matter density. The CRs diffuse into both the disk and the halo. 
The diffusion coefficient is taken to be rigidity dependent and 
position independent such that K{R) - [iK(){RIRQf. The number 
density Nj of the nucleus j is a function of the kinetic energy 
per nucleon, E, and the position (r, z). The steady-state transport 
equation can be written as 



((Wf-Kj'^^)Nj^Sf 



(22) 



The loss term, ^'°', describes the decay rate of unstable nuclei, 
Pf"* = l/ijLTj) (where -yi is the usual Lorentz factor) and the 



total destruction rate for collisions in the disk, 2h6iz)t' . The 
source term, >S'°', is the sum of contributions from SNRs (the 
DSA solution of Eq. 18), and the secondary production in the 
ISM from A;-type progenitors 

Sy ^2h6(z)s(r)Sf\E), (23) 

2li6{z)Y,{ft7+nf)Nk. (24) 

k>j 

The function s{r) expresses the SNR radial distribution in the 
disk, that we assume to be uniform. For primary CRs, 5'^^'' 
is normalized by the Y constants of Eq. 12. In our study, sec- 
ondary nuclei may have a non-zero S'^'^^ term. The term f J!"' = 

{jLTkjY^ describe the contributions k ^> j from unstable pro- 



ST 



genitors of lifetime r^.y. In Eq. 24, P = ;6yc«isnicr™ 

f/^" = PjCiVi^mO'"^ ■ The conditions Nj s at the halo bound- 
aries and the continuity condition across the disk completely 
characterize the solution of Eq. 22. The full solution is reported 
in Maurin et al. (2001). We again solve the transport equations 
for all the CR nuclei following their top-down fragmentation 
sequence, plus a second iteration to account for the "^Be^'^B 
decay. The differential fluxes as a function of kinetic energy per 
nucleon E are obtained from 



dNj Be 



(25) 



where the equilibrium solutions are computed at the Solar 
System position (r, z) - (8.5 kpc, 0). The basic DHM predictions 
can be seen, for illustrative purpose, in the one-dimensional limit 
rmax -^ ^- The solution for a pure primary CR is given by 

cdsa c dsa 

N„ ^ ^ ,. , ~ — , (26) 

^ KIL + r^^ K/L 
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where the spallation rate is neglected for simplicity. The effect 
of the propagation in steepening the spectrum is clearly evi- 
dent: for a source spectrum S'^'^'^iE) oc £ '' and a Galactic dif- 
fusion coefficient of the type K{E) oc E^, the model predicts 
that N{E) oc E^^^^. To resolve the two parameters v and 6, one 
has to consider pure secondary species. The solution for a one- 
progenitor secondary CR is given by Eq. 26, with the replace- 
ment ofS'^''' with 5'^" = fsp^Np, so that the ratio NJNp oc L/K 
allows the simultaneous determination of 6 and Kq/L. These 
simple trends are valid for the reference model, i.e. when only 
primary CRs have a source term. In the case of a secondary SNR 
component, depending on its intensities, the parameter determi- 
nation may be more complicated. 



4. Analysis and results 

We now review the model parameters and test the reference 
model setup. We then analyze the secondary CR production and 
re-acceleration in SNRs in some simple scenarios. 



4.1. Model parameters 

The DSA mechanism of Sect. §2 provides power-law spectra 
^dsa ^ J7-V ^jj.jj ^ unique spectral index v - a - 2 for all the 
primary CRs, where a, in turn, is linked to the compression ratio 
r, which is specified by 6 and the observed log-slope y. To match 
7 * 2.7 with 6 < 0.7, one has to adopt a compression factor of 
r < 4, in contrast to the value r - 4 required for strong shocks. 
Despite this tension between DSA and observations, we regard 
r as an effective quantity describing the compression ratio actu- 
ally felt by the particles, which is not necessary related to the 
physical strength of the SNR shocks (Ptuskin et al., 2010). The 
diffusion coefficient around the shock is taken to be Bohm-like, 



D 



_pc_ 
3ZB- 



The ambient magnetic field B may reach ~ 100/iG 



or more, because of amplification effects, except for the very 
late SNR evolutionary stages, where the magnetic field may be 
damped (B < jjG). In our steady-state description, the SNR pa- 
rameters have to be considered as effective time-averaged quan- 
tities representing a more complex situation where the shock 
structure evolves with time and may be influenced by the back- 
reaction of accelerated CRs. The average shock speed ui is of 
the order of 10* cms"'. The upstream gas density, «i, is poorly 
known and may well vary from ~ 10"^ to ~ 10 cm"^, depending 
on the SNR progenitor star or its local environment. The SNR 
explosion rate per unit volume is expressed as a surface density, 
2h'R,„„ that we fix to 25 Myr"' kpc^^ (Grenier, 2000). 

The parameters describing the interstellar diffusion coeffi- 
cient are fixed to 6 - 0.5 and Kq = 0.089 kpc^ Myr' (see 
Sect §4.2). Below the reference rigidity, Rq = 4 GV, we set 6=0. 
However our analysis is always applied to rigidities R > Rq. The 
halo radius is rmax = 20 kpc and its half-height is L = 5 kpc. 
As per the propagation in the ISM, the quantity that enters the 
model is surface density h x «ism, where we take /i = 0.1 kpc and 
"ism = 1 cm"^. We assume a composition of 90% H + 10% He 
for the ISM gas density, «ism, and that this composition, on av- 
erage, is also found in the SNR background media. In addi- 
tion, we include the solar modulation effect, though it is relevant 
only below a few GeVnucleon '. The modulation is described 
in the force-field approximation (Gleeson & Axford, 1968) by 
means of the parameter 0, taken to be 500 MV, to characterize a 
medium-level modulation strength. Our nuclear chain starts with 
Zmax - 14 and processes all the relevant isotopes down to Z = 3. 
Nuclei with Z > 14 do not contribute significantly to the Li- 



Be-B abundances. The spallation cross-sections are taken from 
Silberberg et al. (1998). The cross-sections on He targets are ob- 
tained by means of the algorithm presented in Ferrando et al. 
(1998). The reference model parameters are listed in Table 1. 



Table 1. Source and transport parameter sets. 



Acceleration parameters 


Propagation parameters 


lit lO'^cms-' 


^0 


0.089 kpc^Myr-' 


B 50.0/^0 


S 


0.50 


r 2.7 


^0 


4GV 


ii[ 1 cm ^ 


L 


5 kpc 


7 14 


''max 


20 kpc 


Tsnr 20 kyr 


h 


0.1 kpc 


??snr 125Myr-'kpc-2 


^sm 


1 cm-3 


Rini 1 GV 





0.5 GV 



4.2. Reference model 

Before analyzing the impact of SNR production and re- 
acceleration of secondary CRs on the parameter determina- 
tion, we test the reference model predictions for the parame- 
ters in Table 1. Predictions at Earth for the CR elemental spec- 
tra Li, Be, B, C, N, O, Ne, Mg, and Si are presented in Fig. 1, 
where the total spectra (solid lines) are shown together their 
secondary component (dashed lines) arising from collisions in 
the ISM. The key quantities for propagation, Kq and S, are 
determined from the B/C ratio above 2 GeVnucleon"', using 
all the data reported in the past two decades, i.e., from the 
space-based experiments HEA03-C2 (Engelmannet al., 1990), 
CRN (Swordy et al., 1990), and AlVIS-Ol (Aguilaret al., 2010), 
and the balloon-borne projects CREAM (Ahn et al., 2009) and 
ATIC-2 (Panov et al., 2007). The primary nuclei spectra were 
normalized using data from CREAM (for C, N, O, Ne, IVIg, and 
Si) and HEA03-C2 (for all elements). Given the Kq - L degen- 
eracy (see Sect. §3), in the following we adopt the quantity K^)/L 
as the physical parameter, where the halo height L is fixed at 
5 kpc. 

As apparent from the figure, the reference model calculations 
give a good description of the CR elemental spectra within the 
precision of the present data. We note that, under this pure dif- 
fusion model, the B/C data between ~ lOGeV and ~ 1 TeV per 
nucleon suggest that S ~ 0.4, while the data at lower ener- 
gies (~ 1-100 GeVnucleon"') favor higher values {6 ~0.6). 
These uncertainties are related on both the model unknowns at 
~ GeV/n energies and the lack of data at > 100 GeVnucleon"' 
(Maurin et al., 2010). We also note that the reference model 
predictions are insensitive to the source parameters ni, u\, and 
B: the source spectra are specified only by the effective com- 
pression ratio r (via y and 6) and the abundance constants Y. 
This setup is equivalent to that of many diffusion models, e.g. 
Maurin et al. (2001), that make use of rigidity power-law param- 
eterizations as source functions. 

4.3. SNR models 

Similarly to Morlino (2011), we consider two ideal situations 
represented by Type I/a (important for fragmentation. Sect. §4.4) 
and core-collapse supemovae (important for re-acceleration. 
Sect. §4.5). In the Type I/a scenario, the supernova explodes 
in the regular ISM with typical density and temperature «i ^ 
1 cm"^ and To = 10"* K. In the core-collapse scenario, the SNR 
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expands into a hot diluted bubble (ni <K 1 cm"^ and To » 10"* K) 
that may be generated by either the progenitor's wind or by pre- 
vious SNR explosions that occurred in the same region. In both 
scenarios, the circumstellar densities are assumed to be homo- 
geneous and constant during the SNR evolution. Two SNR evo- 
lutionary stages are relevant to our study: the ejecta-dominated 
(ED) phase, when the shock front expands freely and accumu- 
lates the swept-up mass in the SNR interior, and the Sedov- 
Taylor (ST) expansion phase, which is driven by the thermal 
pressure of the hot gas. The phase transition ED-ST occurs at the 
time Tst, when the swept-up mass equals the mass of the ejecta 
Mej. The CR acceleration ceases at the time Tsnr- 



Table 2. Case studies of Type I/a and core-collapse SNRs. 



SNR model 



ii (cm ) Ta (yr) 



Ml (cms ') 



type I/a 



I/a#l 
I/a #2 
I/a #3 



0.5 
1.5 
3.0 



330 
230 
183 



1.3x10'* 
1.0x10** 
8.9x10^ 



core CC #1 


0.003 


1829 


3.5x10* 


collapse CC#2 


0.01 


1225 


2.8x10** 


CC#3 


0.1 


568 


1.8x10** 



For all the models, we assume Esm- = 5-10^' erg, /Wgj - 4-Mq, 
and Tsnr = 20kyr, where E^ni is the SNR explosion energy (not 
converted into neutrinos) and Mq is the solar mass. The cases of 
different T^m- values are considered in Sect. §4.6. During the ED 
stage, the SNR radius grows with a constant rate R^hit) - u\t, 

at the speed ui - (2£'snr/Mejj ~ 10^ cms"', until it reaches 
the swept-up radius R^v, = RshiTsi)- For a given SNR model 
characterized by «i and Tsm, we parametrize the shock evo- 
lution (/?sh and Mi) using the self-similar solutions derived in 
Truelove & McKee (1999) for a remnant expanding into a ho- 
mogeneous medium. These solutions connect smoothly the ED 
phase (/?sh oc with the ST stage (R^h <x t^'^)- The CR accel- 
eration ceases at time Tsnr, from which we compute the average 
velocity mi s ^sh(Tsnr)/T"sm- Thus, wc usc M~i as an input parame- 
ter for our steady-state DSA calculations (Sect. §2.1) to compute 
the spectra for all the CR elements. The SNR models considered 
are listed in Table 2. 

We always assume that the total CR flux is produced by 
SNRs of only one type: for each SNR model, the parameters 
employed have to be regarded as effective ones representing the 
average population of CR sources. We note that this simplified 
breakdown is somewhat artificial, because the total CR flux may 
be due to a complex ensemble of contributing SNRs. 

4.4. Secondary CR production in Type I/a SNRs 

The secondary production of CRs is relevant for SNRs that ex- 
pand into ambient densities of the order of «i ~ 1 cm"^, where 
the quantity «i represents the average SNR background density. 
Such a value may be higher than that of the average ISM, owing 
to, e.g., contributions from SNRs located in high density regions 
of the Galactic bulge, inside the dense cores of molecular clouds 
or SNRs expanding into the winds of their progenitors. 

From Eq. 16, one sees that the secondary CR flux emitted 
by SNRs has two components. By analogy with Blasi & Serpico 
(2009) and KacheMeB et al. (2011), these are referred to as ^ 
and S. The yi-term, proportional to /o, describes the particles 
that are produced within a distance ~ D/u of both the sides 
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Fig. 2. Energy spectra of secondary elements Li, Be, and B. The source 
components from fragmentation occurring inside SNRs are split into the 
^-term (dot-dashed lines) and S-term (dotted lines). The dashed lines 
indicate the ISM-induced components. The solid lines represent the to- 
tal spectra. Source parameters are reported in the text. The propagation 
parameters are as in Table 1 . Data as in Fig. 1 . 



of the shock front and are still able to undergo DSA. From 
Eq. 17, the ^-spectrum is / ~ p^"^^, reflecting the spectrum 
of their progenitors (/ ~ p^") and the momentum dependence 
of the diffusion coefficient, D cc p. The S-term, / ~ qo, de- 
scribes those secondary nuclei that, after being produced, are 
simply advected downstream without experiencing further ac- 
celeration. Their spectrum maintains the same behavior as that 
of their progenitors q2 ~ p'" ■ These components are illus- 
trated in Fig. 2 for the spectra at Earth of Li, Be, and B for a 
SNR with «i = 2cm-^ ui ^S-lO^cms^', and B = 0.1//G. 
The figure compares the standard ISM components (solid lines) 
with the source components J?( (short-dashed lines) and S (long- 
dashed lines) within the same propagation parameter set. Both 
the source components are harder than those expected from the 
standard CR production in ISM; in particular, the ^[-term leads 
to increasing secondary-to-primary ratios at high energies (Blasi, 
2009; Mertsch & Sarkar, 2009). While the S-term depends on 
the SNR ambient density «i and its age Tsnr, the ^-term also 
relies on the diffusion properties, as its strength is proportional 
to ~ Y^'^^DIu\. However, the parameter combination ni/(Bu^^) 
should be sufficiently large to ensure an y[-term dominance at 
~ TeV energies, which can be realized only in the latest evolu- 
tionary stage of a SNR characterized by damped magnetic fields 
(B «clyuG) and low shock speeds (ui <10'*cms '). On the 
other hand, the local flux of stable CR nuclei depends on the 
large-scale structure of the galaxy (of some kpc) and reflects the 
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Fig. 4. Fit results for the parameters 6 and Kq/L of our models with fragmentation in type I/a SNRs. Results are shown for E^i„ =2, 5 and 
lOGeVnucleon"' (top to bottom), for the reference model and for the SNR models I/a # 1, #2 and #3 (left to right) of Table 2. The shaded areas 
represent the 1-, 2- and 3-cr contour limits. The markers "x" indicate the best-fit parameters for each configuration; the^-/rf/ ratio is reported in 
each panel. 



contribution of a relatively large population of SNRs and their 
histories (Taillet & Maurin, 2003). Furthermore, from Eq. 9 and 
Eq. 1 1, the yi-term induces an exponential cut-off at momentum 
p'="' given by ;f(p'^"') ~ 1, which is not observed in present data 
of primary or secondary CR spectra. Since our aim is to estimate 
these effects when the considered SNRs produce all the observed 
CR flux, the associated parameters have to be able to accelerate 
all CR nuclei up to, say, /7™™/Z ~ lO^'GV. Thus, from the re- 
quirement that;^^ < 1 for any p up to /9™'' (see Sect. §2.1), the 
J?l-term is always ineffective at the energies we consider and is 
not discussed below. Given the absence of a clear spectral feature 
in the S-term, the spectral deformation induced by interactions 
in SNRs may be difficult to detected ~ TeV energies, because 
it can be easily mimicked by a different choice of 6 and /To/L. 
This is illustrated in Fig. 3, where the total boron spectrum (solid 
line) is plotted showing its standard component arising from ISM 
collisions (dashed line) and the source component coming from 
hadronic interactions in SNRs (dotted line). The SNR model is 
the I/a #3 of Table 2. We note that the carbon flux also contains a 
small amount of secondary fragments (<5%), produced in both 
ISM and SNRs. The B/C ratio is also plotted for the reference 
model (dashed line) under the same propagation parameter set- 
ting, i.e. when hadronic interactions in SNRs are turned off. The 
effect of including secondary production in the sources translates 
into a slight increase at 100 GeV nucleon"', while it reaches a 
factor 2.5 at 1 TeV nucleon ' and one order of magnitude at 10 
TeV nucleon '. 



From the B/C ratio data, we determined the parameters K^jL 
and S for the Type I/a SNR models of Table 2. We performed a 
X^ analysis using our model interfaced with MINUIT. The data 
were fit above a minimal energy of £■,„!„ = 10 GeV nucleon"', 
as a compromise between the diffusion-dominated regime and 
the availability of experimental data. We also repeated the fits 
down to lower £■„,;„ to test the relevance of low energy effects 
to our model. The results are shown in Fig. 4 for £■,„!„ - 2, 5, 
and 10 GeV nucleon"' (from top to bottom), for reference model 
and I/a models of Table 2 (left to right). The shaded areas repre- 
sent the 1-, 2-, and 3-cr contour limits of the x^- We stress that 
these parameter uncertainties are those arising from the fits and 
that they are contextual to our models. Owing to the complex- 
ity of the physics processes involved together with the possible 
lack of knowledge of several astrophysical inputs, the actual pa- 
rameter uncertainties may be much larger (Maurin et al., 2010). 
For instance, the published values of 6 vary widely from ~0.3 
to ~ 0.7. The markers describe the best-fit parameters for each 
configuration. The x^ values reported in each panels are divided 
by the degrees of freedom df - 26, 21, and 16 for the considered 
energy thresholds. It can be seen from Fig. 4 that the source com- 
ponent has a little effect for model I/a #3 («i = 0.5 cm""*). When 
denser media were considered, the secondary source component 
was found to flatten the B/C ratio, so that higher values of 6 
were required to match the data. This trend is clearly apparent in 
Fig. 4 (from left to right). Similar conclusions, though weaker, 
can be drawn for the KqJL parameter ratio. To first approxima- 
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Fig. 3. Top: individual CR spectra of B and C. Solid lines are the model 
predictions for I/a #3 SNR model of Table 2. Model parameters are as 
in Table 1, except for 6 and Kq/L, which are fitted to data. The boron 
SNR component (dotted line) and the ISM component (dashed lines) 
are reported. Bottom: the B/C ratio from the above model (solid line) 
and when fragmentation in SNR is turned off (dashed line). Data are 
from HEA03-C2 (Engelmann et al., 1990), CREAM (Ahn et al., 2009), 
AMS-01 (Aguilar et al., 2010), TRACER (Obermeier et al., 2011), 
ATIC-2 (Panovetal., 2007), CRN (MuUer et al., 1991), Simon et al. 
(1980), Lezniak & Webber (1978), and Orth et al. (1978). 



tion, B/C oc L/Ko, so that the presence of a SNR component of 
boron requires a larger Kq/L ratio to match the data. 

In summary, for the SNR models considered, the fragmenta- 
tion in SNRs affects the parameter S of ~ 5-15% (and Kq/L of 
~ 2-10%), but these models cannot be discriminated by present 
data because of the large uncertainties in the data. This « i-S de- 
generacy is apparent by the x^ /df-values, which are almost in- 
sensitive to the SNR properties. 

4.5. Re-acceleration in core-collapse SNRs 

The amount of re-accelerated CRs depends on the total volume 
occupied by the SNRs (per unit time) and their explosion rate 
(per unit volume). The fraction of re-accelerated CRs to the 
total background CRs can be roughly estimated as N'^'/N^^ ~ 
Vsnr'^sniTesc, whcrc Vsm is the SNR volumc and Tesc is the 
characteristic escape time of CRs in the Galaxy. At a few 
GeV nucleon"', Tesc ~ 2liL/K ~ 5 Myr The Vsm is mainly deter- 
mined by its expansion during the ED phase; the SNR reaches a 
spherical volume Vsw - Mf.j/{mni), where m is the mean mass 
of the ambient gas. Thus N'^/N^^ oc 1/ni, which is an oppo- 
site trend to that of the fragmentation scenario of Sect. §4.4. 
One can see that for a density «i ~ 1 cm""* the re-acceleration 
gives a small contribution to the total CR flux. In contrast. 



Fig. 5. Top: individual CR spectra of B and C. Solid lines are the model 
predictions for CC #2 SNR model of Table 2. Model parameters are as 
in Table 1 except for S and Kq/L, which are fitted to data. The boron 
SNR component (dotted line) and the ISM component (dashed lines) 
are reported. Bottom: the B/C ratio from the above model (solid line) 
and when re-acceleration is turned off (dashed line). Data as in Fig. 3. 



for «i <0.01cm"^, the re-acceleration fraction grows signifi- 
cantly (> few percent). However it is also important the subse- 
quent ST phase, where the SNR shock expands adiabatically as 
^sh(0 0^ t^^^, slowing down at the rate ui{t) - t"^^^ . 

Using the SNR parameters of Table 2, we computed the re- 
accelerated CR spectra as in Sect. §2.2, using Q"^'^'^ - f'°^{p)6{x), 
where f°°{p) - j;^- Since CRs are already supra-thermal, 

we assumed that all CR particles above p'"J are suitable for 
(re-)undergoing DSA. We note that N(p) is the DHM solu- 
tion of Eq. 22 that, in turn, is fed by the total DSA spectra. 
Hence, we solved the DSA and DHM equation systems iter- 
atively. At the first iteration, only the standard injection term 
was considered (Eq. 12) to compute the interstellar flux N for 
all CR nuclei. The subsequent iterations made use of the previ- 
ous DHM solutions, N, to update the terms QP" and Q'^'^'^ and 
to re-compute the total interstellar fluxes. The procedure was 
iterated until the convergence was reached. At each iteration, 
the injection constants, Y, were re-adjusted. The resulting CR 
flux (standard plus re-accelerated) was therefore determined by 
Eq. 18 and is fully specified by the source parameters n\ and 
Tsni. In practice, we found that five iterations ensure a stable so- 
lution. The effect of re-acceleration is shown in Fig. 5 for the 
SNR model CC #2 of Table 2. At energies of ~ 1 TeV nucleon ', 
the re-accelerated component dominates over the ISM-induced 
component for secondary nuclei. It should be noted that the 
sources of re-accelerated CRs may have a complex spatial dis- 
tribution depending on the SNR spatial profile. In our model, we 
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Fig. 6. Fit results for the parameters S and Kg/L of our models with re-acceleration in core-collapse SNRs. Results are shown for iJ„,i„ = 2, 5, and 
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represent the 1-, 2-, and 3-cr contour limits. The markers "x" indicate the best-fit parameters for each configuration; the x^/df ratio is reported in 
each panel. 



used a uniform distribution, s(r) = 1, which does not limit the 
predicting power of diffusion models as long as the key param- 
eters are regarded as effective quantities tuned to agree with the 
data (Maurin et al., 2001). However, further elaborations would 
require more refined descriptions. In Fig. 6, we plot the fit re- 
sults for the parameters 6 and KqJL for the core-collapse SNR 
models of Table 2 and the reference model. Compared to the sce- 
nario of Sect. §4.4, the results are less trivial to interpret because 
the background CR flux that is subjected to re-acceleration it- 
self depends on the parameters S and KqIL. As consequence of 
this non-linearity, the KqJL best-fit values results are less sensi- 
tive to n\. The results for 6 are qualitatively similar to those of 
Sect. §4.4, showing an opposite dependence on «i. For the SNR 
model CC #l{ii\ - 0.003 cm"-', left column), the source compo- 
nent dominates the secondary CR flux at ~ lOOGeVnucleon ', 
so that, at higher energies, the B/C ratio becomes appreciably 
flat. The effect becomes less significant for higher background 
densities, e.g., CC# 1 («i = 0.1 cm""*, right column) where the 
best-fit parameters are close to those arising from the reference 
model fit. As for the scenario of Sect. §4.4, results are limited by 
the sizable uncertainties in the parameters that preclude quan- 
titative conclusions for ^min - lOGeVnucleon"'. Nonetheless, 
the figure shows clear trends, especially for 5. The x^/df val- 
ues reported in each panel indicates that good fits can be made 
for all the considered configurations, though they do not vary 
significantly among the various SNR models. 



4.6. Summary and discussion 

Our breakdown into Type I/a and core-collapse SNR scenarios 
is motivated by the complementary dependence of the two ef- 
fects on «i. As seen in Sect. §4.4 and Sect. §4.5, the CR re- 
acceleration is found to be important for SNRs exploding into 
rarefied media, which are typical of super-bubbles (including 
our own local bubble), while the secondary production in SNRs 
is relevant for ambient densities similar to those of the regular 
ISM. Our calculations of the B/C ratio are in substantial agree- 
ment with the work of Berezhko et al. (2003) in the cases where 
the comparison can be made, though those authors used different 
approaches to model the acceleration as well as the interstellar 
propagation. The fit results for the SNR models of Table 2 and 
the reference model are listed in Table 3. 

Figure 7 summarizes our findings, showing the best-fit pa- 
rameters as functions of the SNR circumstellar density. The 
panel groups (a), (b), and (c) are referred to fits performed at 
different minimal energies, Eatin = 10, 5, and 2GeVnucleon ', 
respectively. For each group, we report 6, Kq/L, and x^/df as 
functions of ni (from bottom to top, solid lines). The two mech- 
anisms are presented separately: the sub-panels on the left-hand 
side show the effect of re-acceleration in CC type SNRs, while 
the right-hand side plots are referred to the secondary produc- 
tion by spallations in Type I/a SNRs. The complementarity of 
the two effects is apparent from the figure. In the region where 
they overlap, «i siO.Scm"^, neither is relevant. The horizon- 
tal (dotted) lines indicate the best-fit parameters for the refer- 
ence model. Their dependence on E,^in resembles that found in 
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Table 3. Summary of the fit results to tfie propagation parameters. See the text in Sect. §4.6. 



Fit 



£ > 2 GeV/n 



£ > 5 GeV/n 



£>10GeV/n 



Model 



Ko/L(kpc/Myr) x^jdf 



A-o/LCkpc/Myr) x^ldf 



/fo/L(kpc/Myr) x^jdf 



Ref. Model 0.50±0.03 0.01781 ±0.00069 15.69/26 0.45±0.04 0.01982 ±0.00149 12.30/21 0.40±0.05 0.02270 ± 0.00270 8.76/16 

SNRI/a#l 0.51 ±0.03 0.01813 ±0.00071 15.27/26 0.46±0.04 0.02008 ±0.00153 12.19/21 0.41 ±0.05 0.02296 ± 0.00280 8.81/16 

SNRI/a#2 0.52 ±0.03 0.01895 ±0.00075 14.34/26 0.48 ±0.04 0.02073 ±0.00165 11.97/21 0.43 ±0.06 0.02357 ± 0.00304 8.97/16 

SNRI/a#3 0.55 ±0.03 0.02022 ± 0.00083 13.34/26 0.52±0.05 0.02164 ±0.00185 11.82/21 0.47±0.06 0.02432 ± 0.00346 9.35/16 

SNRCC#1 0.62 ±0.03 0.01897 ±0.00082 13.45/26 0.60 ±0.04 0.02005 ±0.00179 12.22/21 0.56 ±0.07 0.02217 ±0.00343 10.16/16 

SNRCC#2 0.57±0.03 0.01847 ±0.00076 13.95/26 0.53±0.04 0.02003 ±0.00165 11.97/21 0.49±0.06 0.02258 ± 0.00308 9.33/16 

SNRCC#3 0.52 ±0.03 0.01798 ±0.00070 15.14/26 0.47 ±0.04 0.01988 ±0.00152 12.15/21 0.42 ±0.05 0.02269 ± 0.00279 8.87/16 
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Fig. 7. Best-fit parameters d and Kg/L and the corresponding x'/df as function of «i (solid lines) for models with re-acceleration in core-collapse 
SNRs (CC, left sub-panels) and with hadronic interactions in Type I/a SNRs (I/a, right sub-panels). The panel groups (a), (b) and (c) show the fit 
results for data at £„!;„ = 1, 10, 5, and 2GeVnucleon"', respectively, for T^m = 20 kyr. Panel (d) shows the results for t„„- = 10, 20, 40, 60, and 
80 kyr in the case of Eniin = 2 GeV nucleon"' . The horizontal dotted lines indicate the reference model parameters. 



Di Bernardo et al. (2010), who also considered diffusive reaccel- 
eration models, although we note that we used a different set of 
data for the parameter determination. As discussed, the refer- 
ence model is insensitive to either n\ or other SNR parameters. 
The dashed lines indicate the parameter uncertainties (at one <j 
of CL) arising from the fits. 



It is interesting to note the evolution of the htsi-x^ 
structures when the minimal energy E^in is decreased from 
lOGeVnucleon"' (Fig. 7a) to 2GeVnucleon"' (Fig. 7c). When 
the low-energy B/C data are included in the fits, i\\t x^ Idf dis- 
tribution exhibits two minima. The B/C ratio data at low ener- 
gies favor a slope (6 ~ 0.6), which is somewhat steeper than 
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that observed in high energy data (6 ~ 0.4): these two regimes 
are matched by models of SNRs that emit secondary nuclei. In 
some previous studies, e.g. Trotta et al. (201 1), it has been found 
that the secondary-to-primary ratios can be reproduced well at 
all energies using (5=1/3 and a strong diffusive reacceleration 
(the interstellar Alfvenic speed is of the order of ~ 30kms '). 
However, this cannot be satisfactorily reconciled with the use 
of pure power-law functions for the CR sources. On the other 
hand, the trends we observe suggest a possible role of SNR- 
fragmentation or re-acceleration in reconciling the low-energy 
B/C data with those at higher energies in a pure diffusion sce- 
nario. 

As we have stressed, the physical effects discussed in this 
work should be tested at high energies, where most of the com- 
plexity of the low-energy CR propagation can be neglected. 
Owing to the scarcity of CR data above lOGeVnucleon ', the 
parameter constraints reported in Fig. 7a do not allow us to 
make any firm discrimination among the different SNR models. 
However, the main trends are apparent. On the propagation side, 
all the non-standard scenarios point toward larger values for 6, 
which is in some tensions with the predictions for the interstellar 
turbulence (Strong et al., 2007) and with CR anisotropy studies 
(Ptuskin et al., 2006). On the acceleration side, large values of 6 
reduce the source spectral index closer to the value v -2, which 
is favored by the DSA theory for strong shocks. 

In all these scenarios, the acceleration ceases at Tsnr - 20 kyr, 
which may not be the case given their different SNR evolu- 
tionary properties. For instance, since the ST phase duration 
scales as n^ ' (Truelove & McKee, 1999), one may expect core- 
collapse SNRs to have longer Tsnr than Type I/a SNRs. However, 
the parameter Tsnr represents the time for which the SNR is ac- 
tive as a CR factory and it can be extremely difficult to estimate. 
Thus, in Fig. 7d, we give the fit results for different values of 
Tsnr from lOkyrto 80 kyr (Zsmin =2GeVnucleon"'). The effect 
of using different Tsnr is clear. The longer the time for which the 
SNR is active, the larger the fragments produced in its interior. 
In practice, the secondary CRs production in Type I/a SNRs is 
characterized by the product «iTsnr- For re-acceleration, a longer 
lifetime allows the SNR to occupy larger volumes. To first ap- 
proximation, the intensity of re-accelerated nuclei increases as 
~ Tsnr/ni. As shown in the figure, for longer values of Tsnr, the 
re-acceleration effect also becomes important for relatively high 
density media. 



5. The projected AMS-02 sensitivity 

We switch now to some estimations for the AMS experiment', 
which is devoted to direct measurements of Galactic CRs across 
a wide range of energy. The prime goals of the AMS project are 
the direct search for anti-nuclei and the indirect search for dark 
matter particles. The first version of the experiment, AMS-01, 
operated in a test flight on June 1998. The final version of exper- 
iment, AMS-02, was successfully installed in the International 
Space Station on May 2011 and will be active for at least ten 
years. AMS-02 is able to identify CR elements from Z = 1 to 
Z = 26 and to determine their energy spectra from ~ 0.5 GeV 
to ~ 1 TeV per nucleon with unprecedented accuracy. We esti- 
mate the AMS-02 capabilities in determining the CR propaga- 
tion properties for the considered scenarios. 



5.1. Projected data 

The AMS-02 sensitivity to CR nuclei measurements is studied 
by the generation of mock data for a given input model. The 
number of y'-type particles recorded by AMS-02 at the kinetic 
energies between Ei and £2 is given by 



AA^,= r 

Je, 



'>j(E)-Sj-gj-rj-dE, 



(27) 



http : //www . ams02 . org 



where (pj is the input spectrum, Qj is the detector geometric fac- 
tor, Sj is the detection efficiency, and Tj is the exposure time. All 
of these quantities are in general energy-dependent and particle- 
dependent. The relevant quantity for our estimates is the expo- 
sure factor f = &0T, which we assume to be both energy and 
particle independent. We consider the cases of ;F =12, 36, and 
ISOm^srday. For values of, e.g., @ =0.45m^sr and £ =90%, 
our choices correspond to one month, three months and one 
year of time exposures, respectively. We adopt a log-energy bin- 
ning using nine bins per decade between 10 GeV and 1 TeV per 
nucleon. The AMS-02 mock data for the B and C fluxes and 
their ratios are shown in Fig. 8a and 8b. The CR fluxes, <^b 
and (pc, are calculated using the SNR model CC#2 (Table 2) 
as an input model. This "true" model is characterized by the 
SNR parameters n\ = 0.01 cm -^ and Tsnr =20 kyr, and trans- 
port parameters KqJL = 0.01847 kpcMyr"' and 6 =0.57. From 
Eq. 27, we compute the statistical error in each B/C data point 
as 1/ y/ANs + 1 / yJANc. Nonetheless, CR measurements are also 
affected by systematic errors, which become increasingly impor- 
tant as the precision increases with the collected statistics. 

5.2. Discrimination power 

We fit the B/C ratio mock data leaving A'o/L and 6 as free pa- 
rameters. These parameters are determined within both the ref- 
erence model (re-acceleration off) and the "true" re-acceleration 
model CC#2 of Table 2. As shown in Fig. 8b, both the mod- 
els can be tuned to reproduce the AMS-02 mock data, but they 
exhibit different functional shapes and deviate at high energies. 
The contour plots in panels (c) and (d) correspond to the best- 
fit parameters of the two models. Contour levels are shown for 
3-cr uncertainty levels corresponding to the three exposure fac- 
tors ;F. As expected, the re-acceleration model fit (d) returns 
the correct parameters, while the reference model fit (c) mises- 
timates the parameters because of inaccurate assumptions about 
the source properties. In fact, when the reference model is forced 
to describe the data, the spectral distortion induced by the re- 
acceleration is mimicked by the use of a lower value for 6. Given 
the precision of the AMS-02 data, this represents the dominant 
"error" in the parameter determination. As apparent from the fig- 
ure, the AMS-02 data place tight constraints on the propagation 
parameters. For instance, 6 is determined to a precision better 
than < 10% within a 3-cr uncertainty level. The 6-ni degeneracy 
may be lifted as in Castellina & Donato (2005), i.e., by a statis- 
tical test to discriminate between the two fits. As long as only 
statistical errors are considered, we find that the discrimination 
between the two scenarios is always possible for the three con- 
sidered exposures at 90 % of CL. The effect of systematic errors 
in the data is shown in Fig. 8e, where we plot the AMS-02 dis- 
crimination probability versus the relative systematic error Our 
calculation assumes constant systematic errors (added in quadra- 
ture to the statistical ones), but these considerations also hold for 
energy-dependent systematic errors if their energy rise is less 
pronounced than the statistical errors. The solid, dashed, and 
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Fig. 8. AMS-02 mock data for the elemental fluxes ipg and ipc (a) and their ratio (b), using the input model CC#2 of Table 2 and assuming a 
detector exposure factor F= 150 m^ srday. The error bars are only statistics. The constraints to the transport parameters provided by the AMS-02 
mock data are reported by the 3-it contour levels for exposure factors of F = 12, 36, and 150 m^ srday in (c) and (d). The data are fitted within the 
reference model (star in panel c and dotted line in panel b) and within the model CC#2 (cross in panel d and solid line in panel b). The AMS-02 
discrimination probability between the two models as a function of the systematic error in the measurement is shown in (d) for F= 12, 36, and 
150m- srday. The systematic errors are assumed to be energy-independent. 



dotted lines represent the cases of 'F = 12, 36, and 150 m^ srday, 
respectively. To achieve a discrimination of 90% CL, the system- 
atic error has to be smaller than ~ 4%, 8%, and 10% for the three 
considered exposures. A 95% CL requirement also needs 'F to 
be larger than 12m- srday. We consider these requirements as 
reasonable for AI\/lS-02, because the measurements of elemental 
ratios are only mildly sensitive to systematic errors. 

Similar conclusions can be drawn for models with fragmen- 
tation in SNRs. In this case, we have explored a large region of 
the parameter space ni-6, with 6 =0.3-0.8 and «i =0-5cm"^. 
Our estimate was carried out as follows. For each {ni,S} param- 
eter combination, we determined K^)/L from fits to the existing 
B/C ratio data. Then we defined the true model using {ni, Tsm) as 
source parameters and {6,Ko/L} as transport parameters. From 
the true model, we generated the AMS-02 mock data for a given 
exposure factor, ;F, and a 5 % systematic error Thus, we re- 
fit the mock B/C ratio, leaving Kq/L and 6 as free parameters, 
within both the reference model and the true SNR scenario (frag- 
mentation specified by «] ). Finally, we estimate the AI\/lS-02 dis- 
crimination probability for the two models. The shaded areas of 
Fig. 9 indicate the parameter region where the AMS-02 discrim- 
ination succeeds at 95% CL for T - 12, 36, and 150m^ srday. 
The figure shows that AMS-02 is sensitive to a large region of 
the parameter space, except for small «i values (small secondary 
SNR component) and/or small S values (hard ISM component), 
when the intensity of the secondary source component is too 
weak to induce appreciable biases in the propagation parameters. 
This is also the case for Kolmogorov-like diffusion (S - 1/3), 
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Fig. 9. AMS-02 discrimination power for models with secondary pro- 
duction in SNRs. Each point in the (iii,S)-p\ane represents an input 
model with fragmentation inside SNRs with Tj,,,. = 20 kyr. The parame- 
ter Ko/L is taken to match the existing B/C ratio data. The shaded areas 
cover the parameter region where AMS-02 is sensitive at 95% CL for 
the exposure !F = 12, 36, and 150 m^ srday. The systematic eiTors are 
assumed to be 5% of the measured B/C and constant in energy. 
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which, however, is disfavored by our analysis of the real data. 
These considerations can be much strengthened if one consid- 
ers the independent constraints that may be brought by other 
AMS-02 data such as, for example, the ratios p/p, Li/C, F/Ne, 
or Ti/Fe. In summary, our estimates show that AI\/lS-02 performs 
wel in determining the CR transport parameters, providing tight 
constraints and considerable progress in understanding the CR 
acceleration and propagation processes. 

6. Conclusions 

We have studied the CR propagation physics under the scenar- 
ios where secondary nuclei can be produced or re-accelerated 
by Galactic sources. We have considered the processes of sec- 
ondary productions inside SNRs and re-acceleration of back- 
ground CRs in strong shocks. The two mechanisms complement 
to each other and depend on the properties of the local ISM 
around the expanding remnants. The secondary production in 
SNRs is significant for dense background media, «i > 1 cm"^, 
while the amount of re-accelerated CRs is relevant to SNRs ex- 
panding into rarefied media, «i < 0.1 cm"^. The consequence 
of both mechanisms is a slight flattening of the secondary-to- 
primary ratios at energies above ~ lOOGeVnucleon"'. For the 
B/C ratio, the increase may be a factor of a few at 1 TeV 
nucleon ' and reach an order of magnitude at lOTeV nucleon '. 
Modeling these effects introduces an additional degeneracy be- 
tween the source and the transport parameters. The diff'usion co- 
efficient index S determined from the B/C ratio measurements 
above ~ 10 GeV nucleon ', was found to be underestimated by 
a factor of > 15 % if the underlying model did not account for 
the hadronic production in SNRs with « i > 2 cm"^ or for re- 
acceleration with «i <0.02cm"^. Nonetheless, the current un- 
certainty in S is much larger as the existing data suffer for a lack 
of precision at £ > 10 GeV nucleon '. We have shown that this 
degeneracy may be at least partially broken with data collected 
by high precision experiments such as AMS-02. Were propaga- 
tion in the Galaxy to be described by a Kolmogorov spectrum 
(5=0.33), it would not be misunderstood with possible source 
effects described in this work, because these are expected to pro- 
duce small distortions in the hard B/C ratio. On the other hand, 
we have shown that for 6 ~ 0.4 - 0.8 an AMS-02 like experi- 
ment will be able to discriminate pure propagation trends from 
a source contribution. Data around TeV nucleon ' energies will 
be clue at this aim. Systematic errors that can be contained to 
the ~ 10% level will not prevent a clear discrimination between 
the reference model and the scenarios with SNR components of 
secondary CRs. 

Data from single elements and antiprotons will help to iden- 
tify the possible effects studied in this research. Closer in- 
spections, including the revision of the role of convection and 
diffiisive reacceleration, will require more data at high ener- 
gies, which may be released soon by a number of ongoing 
experiments. The long-duration balloon projects CREAM and 
TRACER, and the space missions AMS-02 and PAMELA are 
currently operating with unprecedented sensitivities and energy 
ranges. Their data will provide valuable pieces of information 
about CR acceleration and propagation physics. 
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